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Baryon and quark superfluidity in the cooling of neutron 
stars are investigated. Observations could constrain combi- 
nations of the neutron or A— hyperon pairing gaps and the 
star's mass. However, in a hybrid star with a mixed phase of 
hadrons and quarks, quark gaps larger than a few tenths of 
an MeV render quark matter virtually invisible for cooling. If 
the quark gap is smaller, quark superfluidity could be impor- 
tant, but its effects will be nearly impossible to distinguish 
from those of other baryonic constituents. 

PACS number(s): 97.60. Jd. 21.65.+f, 95.30.q 

Pairing is unavoidable in a degenerate Fermi liquid if 
there is an attractive interaction in any channel. The 
resulting superfluidity, and in the case of charged par- 
ticles, superconductivity, in neutron star interiors has a 
major effect on the star's thermal evolution through sup- 
pressions of neutrino (u) emission processes and specific 
heats 0J^]. Neutron (n), proton (p) and A-hyperon su- 
perfluidity in the x 5o channel and n superfluidity in the 
3 P2 channel have been shown to occur with gaps of a few 
MeV or less J^,||. However, the density ranges in which 
gaps occur remain uncertain. At large baryon densities 
for which perturbative QCD applies, pairing gaps for like 
quarks have been estimated to be a few MeV ||]. How- 
ever, the pairing gaps of unlike quarks (ud, us, and ds) 
have been suggested to be several tens to hundreds of 
MeV through non-perturbative studies Q kindling inter- 
est in quark superfluidity and superconductivity |?],|| and 
their effects on neutron stars. 

The cooling of a young (age < 10 5 yr) neutron star is 
mainly governed by v— emission processes and the specific 
heat Q . Due to the extremely high thermal conductivity 
of electrons, a neutron star becomes nearly isothermal 
within a time t w 1 — 100 years after its birth, depending 
upon the thickness of the crust @]. After this time its 
thermal evolution is controlled by energy balance: 

ig.o r §-L,-L. + n, (i) 

where E t h is the total thermal energy and CV is the spe- 
cific heat. L 7 and L v are the total luminosities of photons 
from the hot surface and vs from the interior, respec- 
tively. Possible internal heating sources, due, for exam- 
ple, to the decay of the magnetic field or friction from 
differential rotation, are included in H. Our cooling sim- 
ulations were performed by solving the heat transport 



and hydrostatic equations including general relativistic 
effects (see [|J). The surface's effective temperature T e 
is much lower than the internal temperature T because 
of a strong temperature gradient in the envelope. Above 
the envelope lies the atmosphere where the emerging flux 
is shaped into the observed spectrum from which T e can 
be deduced. As a rule of thumb T e /10 6 K w ^T/10 8 K, 
but modifications due to magnetic fields and chemical 
composition may occur. 

The simplest possible v emitting processes are the di- 
rect Urea processes f\ + I — *■ f% + i>t , f% — ► /i + 1 + VI, 
where f\ and ji are either baryons or quarks and £ is 
either an electron or a muon. These processes can oc- 
cur whenever momentum conservation is satisfied among 
/i , /2 and I (within minutes of birth, the v chemical po- 
tential vanishes). If the unsuppressed direct Urea process 
for any component occurs, a neutron star will rapidly 
cool because of enhanced emission: the star's interior 
temperature T will drop below 10 9 K in minutes and 
reach 10 7 K in about a hundred years. T e will hence 
drop to less than 300,000 K after the crustal diffusion 
time t w Jlp|jrc|] . This is the so-called rapid cooling 
paradigm. If no direct Urea processes are allowed, or they 
are all suppressed, cooling instead proceeds through the 
significantly less rapid modified Urea process in which 
an additional fermion enables momentum conservation. 
This situation could occur if no hyperons are present, 
or the nuclear symmetry energy has a weak density de- 
pendence HQ. The v emisssion rates for the nucleon, 
hyperon, and quark Urea and modified Urea processes 
can be found in fl3[| . 

The effect of the pairing gaps on the emissivities and 
specific heats for massive baryons are investigated in fji|] 
and are here generalized to the case of quarks. The prin- 
cipal effects are severe suppressions of both the emissivity 
and specific heat when T << A, where A is the pairing 
gap. In a system in which several superfluid species exist 
the most relevant gap for these suppressions is the small- 
est one. The specific heat suppression is never complete, 
however, because leptons remain unpaired. Below the 
critical temperature T c , pairs may recombine, resulting 
in the emission of vv pairs with a rate that exceeds the 
modified Urea rate below 10 10 K Jl5| |; these processes are 
included in our calculations. 

The baryon and quark pairing gaps we adopt are shown 
in Fig. Note that gaps are functions of Fermi momenta 
(j>f(i), i denoting the species) which translates into a 
density dependence. For pp(n,p) < 200 — 300 MeV/c, 
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nucleons pair in the state, but these momenta cor- 
respond to densities too low for enhanced v emission in- 
volving nucleons to occur. At higher p_f's, baryons pair 
in higher partial waves. The n 3 J ) 2 gap has been calcu- 
lated for the Argonne Vis, CD-Bonn and Nijmegen I & 
II interactions This gap is crucial since it extends to 
large j?p(n) and can reasonably be expected to occur at 
the centers of neutron stars. For pF{n) > 350 MeV/c, 
gaps are largely uncertain because of both experimental 
and theoretical uncertainties Q. The curves [a], [b] and 
[c] in Fig. |] reflect the range of uncertainty. The p 3 P2 
gap is too small to be of interest. Gaps for the 1 So pair- 
ing of A, taken from Q and shown as dotted curves, are 
highly relevant since As participate in direct Urea emis- 
sion as soon as they appear |l2]| . Experimental informa- 
tion beyond the 1 So channel for A is not available. As for 
£— hyperons remain largely unexplored. The quark (q) 
gaps are taken to be Gaussians centered at pf(q) = 400 
MeV/c with widths of 200 MeV/c and heights of 100 MeV 
[model D], 10 MeV [C], 1 MeV [B] and 0.1 MeV [A], re- 
spectively. The reason for considering quark gaps much 
smaller than suggested in [||^] is associated with the mul- 
ticomponent nature of charge-neutral, beta-equilibrated, 
neutron star matter as will become clear shortly. 
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FIG. 1. Pairing gaps adopted for neutron 1 So and 3 P2, pro- 
ton 1 So, A 1 So, and quarks. The n 3 P2 gaps are anisotropic; 
plotted values are angle- averaged. The A gaps correspond, in 
order of increasing A, to background densities ns = 0.48, 0.64 
and 0.8 fm -3 , respectively. The s-wave quark gaps are 
schematic; see text for details. 



We consider four generic compositions: charge-neutral, 
beta equilibrated matter containing nucleons only (np), 
nucleons with quark matter (npQ), nucleons and hyper- 
ons (npH), and nucleons, hyperons and quarks (npHQ). 
In the cases involving quarks, a mixed phase of baryons 
and quarks is constructed by satisfying Gibbs' phase 
rules for mechanical, chemical and thermal equilibrium 
[Oj. The phase of pure quark matter exists only for 



very large baryon densities, and rarely occurs in our neu- 
tron star models. Baryonic matter is calculated using a 
field-theoretic model at the mean field level |l7| ]; quark 
matter is calculated using either a bag-like model or the 
Nambu-Jona-Lasinio quark model [ pTsj ] . The equation of 
state (EOS) is little affected by the pairing phenomenon, 
since the energy density gained is negligible compared to 
the ground state energy densities without pairing. 

Additional particles, such as quarks or hyperons, have 
the effect of softening the EOS and increasing the cen- 
tral densities of stars relative to the np case. For the 
npQ model studied, a mixed phase appears at the den- 
sity ns — 0.48 fm~ 3 . Although the volume fraction of 
quarks is initially zero, the quarks themselves have a sig- 
nificant pf{q) when the phase appears. The pfs of the 
three quark flavors become the same at extremely high 
density, but for the densities of interest they are differ- 
ent due to the presence of negatively charged leptons. In 
particular, pf(s) is much smaller than pf{u) and pf{<1) 
due to the larger s-quark mass. Use of the Nambu- 
Jona-Lasinio model, in which quarks acquire density- 
dependent masses resembling those of constituent quarks, 
exaggerates the reduction of pp(s). This has dramatic 
consequences since the pairing phenomenon operates at 
its maximum strength when the Fermi momenta are ex- 
actly equal; even small asymmetries cause pairing gaps 
to be severely reduced ||,|ll|. In addition, one may also 
expect p-wave superfluidity, to date unexplored, which 
may yield gaps smaller than that for the s-wave. We 
therefore investigate pairing gaps that are much smaller 
than those reported for the case of s-wave superfluidity 
and equal quark p^'s. 

The introduction of hyperons does not change these 
generic trends. In the case npH, the appearance of hy- 
perons changes the lepton and nucleon pf's similarly to 
the appearance of quarks although with less magnitude. 
While the appearance of quarks is delayed by the exis- 
tence of hyperons, at high densities the j»f's of nucle- 
ons and quarks remain similar to those of the npQ case. 
The existence of either hyperons or quarks, however, does 
allow the possibility of additional direct Urea processes 
involving themselves as well as those involving nucleons 
by decreasing pj?(n) — pf(p)- For the npQ and npHQ 
models studied, the maximum masses are = 1.5M Q , the 
central baryon densities are = 1.35 fm~ 3 , and the volume 
fractions of quarks at the center are = 0.4. 

Cooling simulations of stars without hyperons and with 
hyperons are compared, in Figs. ^ and |^, respectively, to 
available observations of thermal emissions from isolated 
neutron stars. Sources for the observational data can be 
found in pjl ]. However, at the present time, the inferred 
temperatures must be considered as upper limits because 
the total flux is contaminated, and in some cases dom- 
inated, by the pulsar's magnetospheric emission and/or 
the emission of a surrounding synchrotron nebula. Fur- 
thermore, the neutron star surface may be reheated by 
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magnetospheric high energy photons and particles; late- 
time accretion for non-pulsing neutron stars is also pos- 
sible. Other uncertainties arise in the temperature es- 
timates due to the unknown chemical composition and 
magnetic field strength in the surface layers, and in the 
age, which is based upon the observed spin-down time. 
In these figures, the bolder the data symbol the better 
the data. 
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FIG. 2. Cooling of 1.4Mq stars with np matter (continuous 
curves) and npQ matter (dashed and dotted curves). The 
curves labelled as [a], [b], and [c] correspond to n 3 P2 gaps 
as in Fig. [j]; [z] corresponds to zero n gap. Models labelled 
[A] , [B] , [C] and [D] correspond to quark gaps as in Fig. 0; [Z] 
corresponds to zero quark gap. 



The np case is considered in Fig. g, in which solid lines 
indicate the temperature evolution of a 1.4 M Q star for 
quarkless matter: case [z] is for no nucleon pairing at all, 
and cases [a], [b] and [c] correspond to increasing values 
for the neutron 3 P2 gap, according to Fig. |[ The field- 
theoretical model employed for the nucleon interactions 
allows the direct nucleon Urea process, which dominates 
the cooling. The unimpeded direct Urea process carries 
the temperature to values well below the inferred data. 
Pairing suppresses the cooling for T < T c , where T is 
the interior temperature, so T e increases with increasing 
A. If the direct Urea process is not allowed, the range 
of predicted temperatures is relatively narrow due to the 
low emissivity of the modified Urea process. We show 
an example of such cooling (curve [Mc]) using the n 3 P2 
gap [c] for a 1.4M with an EOS Q for which the direct 
Urea cooling is not allowed. 

The other curves in the figure illustrate the effects of 
quarks upon the cooling. The dotted curves [Z] are for 
vanishingly small quark gaps; the dashed curves ([Ah [B], 
[C] and [D]) are for quark gaps as proposed in Fig. [j]. For 



nonexistent ([z]) or small ([a]) nucleon gaps, the quark 
Urea process is irrelevant and the dependence on the ex- 
istence or the size of the quark gaps is very small. How- 
ever, for large nucleon gaps ([b] and [c]), the quark direct 
Urea process quickly dominates the cooling as the nu- 
cleon direct Urea process is quenched. It is clear that for 
quark gaps of order 1 MeV or greater ([B], [C] or [D]) 
the effect of quarks is again very small. There is at most 
a slight increase in the stars temperatures at ages be- 
tween 10 1 to 10 5-6 years due to the reduction of pf(ti) 
and the consequent slightly larger gap (Fig. Q). Even if 
the quark gap is quite small ([A]), quarks have an effect 
only if the nucleon gap is very large ([b] or [c]), i.e., sig- 
nificantly larger than the quark gap: the nucleon direct 
Urea process is suppressed at high temperatures and the 
quark direct Urea process has a chance to contribute to 
the cooling. We find that the effects of changing the stel- 
lar mass M are similar to those produced by varying the 
baryon gap, so that only combinations of M and A might 
be constrained by observation. 
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FIG. 3. Cooling of stars with npH (continuous lines) and 
npHQ matter (dotted lines) for various stellar masses (in 
Mq ). n 3 P2 gaps are from case [clwhile quark gaps, when 
present, are from model [C] of Fig. |l[ 



The thermal evolution of stars containing hyperons has 
been discussed in [^2|,^3| , but we obtain qualitatively dif- 
ferent results here. Hyperons open new direct Urea chan- 
nels: A — > p + e + V e and A + e — > E~ + v e if S _ 's are 
present, with their inverse processes. Previous results 
showed that the cooling is naturally controlled by the 
smaller of the A or n gap. However, this is significantly 
modified when the A gap is more accurately treated. At 
the A appearance threshold, the gap must vanish since 
Pi? (A) is vanishingly small. We find that a very thin 
layer, only a few meters thick, of unpaired or weakly 
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paired A's is sufficient to control the cooling. This effect 
was overlooked in previous works perhaps because they 
lacked adequate zonal resolution. 

In Fig. |3| we compare the evolution of stars of different 
masses made of either npH or npHQ matter. We find 
that all stars, except the most massive npH ones, follow 
two distinctive trajectories depending on whether their 
central density is below or above the threshold for A ap- 
pearance (= 0.54 fin -3 in our model EOS, the threshold 
star mass being 1.28M Q ). In the case of npH matter, 
stars with M > 1.50M© are dense enough so that the A 
gap vanishes and hence undergo fast cooling, while 
stars made of npHQ matter do not attain such high den- 
sities. The temperatures of npH stars with masses be- 
tween 1.3 and 1.5 Mq are below the ones obtained in 
the models of Fig. || with the same n 3 P2 gap [b], which 
confirms that the cooling is dominated by the very thin 
layer of unpaired A's (the slopes of these cooling curves 
are typical of direct Urea processes). Only if the n 3 P2 
gap <C, 0.3 MeV do the cooling curves fall below what is 
shown in Fig. ^. Notice, moreover, that in the mass range 
1.3 - 1.48 M Q the cooling curves are practically indistin- 
guishible from those with unpaired quark matter shown 
in Fig. |^. In these models with npH or npHQ matter, 
there is almost no freedom to "fine-tune" the size of the 
gaps to attain a given T e : stars with A's will all follow 
the same cooling trajectory, determined by the existence 
of a layer of unpaired or weakly paired A's, as long as the 
n 3 P2 gap is not smaller. It is, in some sense, the same 
result as in the case of np and npQ matter: the smallest 
gap controls the cooling and now the control depends on 
how fast the A 1 So gap increases with increasing pp(A). 

Our results indicate that observations could con- 
strain combinations of the smaller of the neutron and 
A— hyperon pairing gaps and the star's mass. Deduc- 
ing the sizes of quark gaps from observations of neutron 
star cooling will be extremely difficult. Large quark gaps 
render quark matter practically invisible, while vanishing 
quark gaps lead to cooling behaviors which are nearly in- 
distinguishable from those of stars containing nucleons or 
hyperons. Moreover, it also appears that cooling obser- 
vations by themselves will not provide definitive evidence 
for the existence of quark matter itself. 
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